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SummMarRY: The solution to the general problem of transferring a rocket 
between two terminals in space with minimum fuel expenditure is explained 
and the results obtained when application is made to a number of particular 
problems of space navigation are described. The mathematical techniques 
which may usefully be employed in the calculation of optimum rocket 
trajectories are exemplified by a method of solving the problem of obtaining 
maximum range from a rocket missile over the Earth’s surface. 


1. Introduction 


In the seventeenth and eighteenth centuries the problems of sea navigation 
acted as a stimulus and a challenge to the scientists of those days and led to many 
enquiries of a mathematical character being undertaken. The results of these so 
fertilised the soil of the whole terrain of mathematics that rich harvests were soon 
gathered from all its fields. The problems of the space navigator are of a far greater 
complexity and, although the mathematical techniques available for solving them 
are more powerful than the primitive methods used by these earlier investigators, it 
is to be expected that the conquest of space and the resulting large demand for 
methods of computing rocket trajectories will have an equally beneficial effect and 
will be productive of many new mathematical ideas. 


There are three main problems of space navigation. Firstly, there is the 
problem of computing the most convenient trajectory to be followed by a space 
ship undertaking a particular interplanetary journey. Among the various possible 
trajectories satisfying the conditions imposed by such factors as the maximum motor 
thrust available, the maximum allowable time for the journey, the minimum 
distance of approach to the sun, and so on, the most satisfactory trajectory will 
clearly be that requiring the least expenditure of fuel. In the early stages, it will be 
necessary to relax all conditions to the greatest degree possible in order to achieve 
fuel economy, while looking forward hopefully to the day when atomic drives 
become a reality and the necessity for fuel economy becomes less pressing. In 
such circumstances, the calculation of trajectories of least fuel expenditure will 
have to be carried out under a variety of conditions imposed by such considerations 
as time schedules and the safety of passengers. For example, it may be a require- 
ment that certain orbits be avoided because of their passing through regions of high 
meteor or cosmic ray intensity. 
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The second main problem is that of computing the position and velocity of a 
rocket in space from observations of celestial bodies and the subsequent derivation 
of the elements of its orbit. This is a survey resection problem in three dimensions 
and, although it is more complex than the corresponding two-dimensional one, the 
existing apparatus of spherical trigonometry will be quite capable of providing a 
solution. The directions of the various bodies of the solar system can be found by 
observing them against the background of the fixed stars and then, knowing the 
co-ordinates of these bodies, those of the rocket are easily derived. Thirdly, there 
is the problem of providing the future navigator with sets of tables, by the aid of 
which he will be able to compute the appropriate manoeuvre necessary to bring his 
ship back on to a pre-calculated track when a divergence has been observed. 


Of these three problems, investigation of the first seems most likely to require 
the development of new mathematical techniques and it is this problem which is 
considered further in the next Section. 


NOTATION 
a= P/p 
A _ see equation (10) 
b see equation (33) 
c rocket exhaust velocity 
see equation (34) 
f(t) acceleration due to motor thrust at time t 
F,, F, components of retardation due to air resistance and gravity 
g acceleration due to gravity 
k see equation (22) 
l,, 1, | direction cosines of direction of motor thrust 
M _ rocket mass 
M, _ value of M on launching 
Pi = xi G@=1, 2) 
P motor thrust 
Gi» Gz rocket’s co-ordinates at “ all-burnt ” 
G1»G2 ~—~rocket’s velocity components at “ all-burnt ” 
R __ range along the Equator 
t time after launching 
T value of ¢ at “ all-burnt” 
X,, X, co-ordinates of the rocket at time t 
X(t) see equation (37) 
yi, % See equations (9) 
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Fig. 1. 
ai, 8, defined after equation (11) 
Kronecker deltas 
6 — inclination of the rocket’s velocity to the horizontal 
¢ inclination of the rocket thrust to the horizontal 
# rate of consumption of propellant 
t= M,/p 


2. Rocket Manoeuvres of Least Fuel Expenditure 


Suppose a rocket is to be navigated between two terminals A and B, its velocity 
at these points being specified. If the motion takes place in vacuo and the gravita- 
tional field is of a general nature, it may be shown"? that the trajectory of minimum 
fuel expenditure comprises a number of null-thrust arcs along which the rocket falls 
freely under gravity, impulsive thrusts from the motor being applied at the junctions 
between these arcs to effect transfer from one to the next (Fig. 1). 


If the motion takes place in the plane of rectangular axes Ox, Oy and the x 
and y components of the gravitational attraction are —f, —g, respectively, two 
components u, v, of a vector quantity called the primer are determined by means 
of the equations 
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Fig. 2. 


The quantities u, v, vary with ¢ along any track of the type specified. The track 
is optimal with respect to fuel expenditure if the quantities u, v are such that 
u?+v? <1 at all points on the track and satisfy the following conditions at each 
junction : — 


(i) u, v are continuous and represent the direction cosines of the direction 
of thrust, 


(ii) wu, v are continuous, 

(iii) “uit +vv=0. 

Corresponding results which take into account aerodynamic forces will be 
found in Ref. 2. 


Particular problems to which this general result may be applied are:— 


(a) The basic manoeuvre of transferring a rocket from one orbit about a 
central attracting body into another not necessarily coplanar with the first. The 
time of transit, the point of departure and the point of arrival, may or may 
not be specified. 
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(a) (0) 


Fig. 3. 


The solution to the simplest case of this problem, namely, transfer between 
coplanar circular orbits, was discovered by Hohmann in 1925"). Hohmann found 
that minimum fuel expenditure is involved if transfer takes place along an elliptical 
orbit tangential to both circular orbits (Fig. 2), impulsive thrusts being applied by 
the motors at each point of contact, A and B, to effect transference into and out of 
the connecting orbit. Although the requirement of impulsive thrusts cannot be 
achieved in practice, there will usually be no difficulty in attaining a very close 
approximation. If the two orbits are those of the Earth and Mars for example, it 
may be calculated that the time which would be spent in the orbit of transfer would 
be about 260 days, as compared with two periods of acceleration of a few minutes 
each. However, account might have to be taken of the fact that the motor thrust 
is limited if transfer between two orbits about the Earth is being considered, for 
then the time of transfer will be comparable with the time of thrust, or if it is 
proposed to use a micro-thrust motor employing a jet of ionised particles. 


The present position in relation to the general problem is that the mode of 
optimal transfer between any pair of coplanar elliptical orbits is known“:* but the 
general three-dimensional manoeuvre remains to be treated. 


(b) It seems probable that a space ship embarking upon an interplanetary 
journey will depart from a circular orbit about the Earth and terminate its passage 
by entering into another circular orbit about the target planet. The problem is 
therefore presented of optimising these manoeuvres for a given pair of planets and 
circular orbits about them. Again, the solution is available“ in the two-dimensional, 
but not in the three-dimensional, case. 


(c) The solution to the previous problem being known, there exists the sub- 


Sidiary problem of determining the most economical manoeuvre whereby a rocket» 
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may escape from a given circular orbit, arriving at infinity with the velocity 
necessary to set it into the orbit of transfer. Assuming the thrust unlimited in 
magnitude, there are two possibilities. If the velocity required at infinity is less 
than a certain critical value, escape should be achieved by a single impulsive thrust 
tangential to the circular orbit (Fig. 3a) causing the rocket to move into a hyperbolic 
orbit. If, however, the velocity at infinity is to exceed this critical value, a more 
complex manoeuvre saves fuel. An impulsive thrust opposing the motion is first 
applied to set the rocket into an elliptical orbit approaching the centre of attraction 
(Fig. 3b). At perigee a second impulsive thrust transfers the ship into a hyperbolic 
orbit along which it proceeds to infinity. This manoeuvre, which was first proposed 
by Oberth"”, becomes more economical the closer the approach that can be made 
to the centre of attraction. By reversing the order of events in these manoeuvres, 
economical entries into circular orbits may be made with any velocity of approach. 
Full details will be found in Refs. 8 and 9. 


If, however, the thrust is limited in magnitude, so that impulsive thrusts cannot 
be approximated, this problem does not admit of so simple a solution. Equations 
are available from which the optimal trajectory may be deduced", but they can- 
not be integrated analytically. Their numerical integration has not yet been 
performed, but it appears that they require the thrust to be offset from the 
direction of motion towards the centre of attraction. 


This section is concluded by describing two other problems relating to optimal 
rocket trajectories. 


When a rocket is perturbed by a planetary body, a transfer of energy takes 
place between it and the body which may be to the advantage of the rocket. Thus, 
by passing close to the Moon, a space ship outward bound from an orbit about the 
Earth to Mars can acquire a considerable amount of energy without any expenditure 
of fuel. The best way of utilising such perturbation effects is not known, although 
there exist a few purely numerical studies, some of which will be found in Ref. 11. 


Finally, there is the most urgent of all astronautical problems, that of calculating 
the optimal trajectory from a ground launching station into a satellite orbit about 
the Earth. This problem is complicated by the existence of aerodynamic forces 
and by the necessity for making allowance for the limited thrust which can be 
provided by the rocket motors. The work of Tsien and Evans?) represents an 
advance along the path leading towards a complete solution to this problem. They 
have calculated the mode of programming the fuel expenditure of a vertically 
ascending rocket necessary to achieve maximum height. The motor thrust was 
assumed unrestricted in magnitude. Clearly a compromise has to be effected 
between a rapid rate of fuel expenditure, leading to high velocities in the lower levels 
of the atmosphere and hence to a fuel wastage due to the excessive work which must 
be done against air resistance, and a slow rate of fuel expenditure, which is 
uneconomical for other reasons. Tsien and Evans’s solution requires that an initial 
impulsive thrust is to be followed by a period of variable finite thrust, at the end 
of which the motor is to be shut down and the vehicle is then to coast to a standstill. 
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They proposed that the initial impulse should be provided by a booster which 
would fall away shortly after the rocket left its launching tower. If c is the jet 
velocity, assuming only that the air resistance is proportional to the square of 
the velocity (i.e. both this force and gravity may vary in any manner with the 
height), optimal programming has been achieved if the speed of the rocket during 
the phases when the motor is operating is Ac, where 


air resistance at jet speed 


The booster is supposed to accelerate the rocket initially to a speed in conformity 
with this equation. A small rocket can be programmed to move in the manner 
required by this equation, but since the rocket weight increases as the cube of the 
linear dimensions, whereas the air resistance only increases as the square, the value 
of A determined by equation (2) increases with the rocket size and there proves to be 
no possibility of a really large rocket following this optimal schedule. It is suggested 
that to reach maximum height with such a rocket, as close an approximation as 
possible should be made to the optimal solution, by first boosting to the maximum 
degree possible and thereafter operating at maximum thrust. For the similar three- 
dimensional problem of placing a large rocket in a satellite orbit, it is to be expected 
that the optimal solution will require a thrust beyond that attainable. Taking into 
account the limited thrust available, optimal conditions will be best approached by 
first boosting to the maximum degree possible and thereafter operating at maximum 
thrust. When just sufficient momentum has been acquired to carry the vehicle to 
the level of the orbit, the motors should be shut down. The rocket then coasts 
to the orbit and a final short burst from the motors will be necessary to effect 
transfer from the elliptical coasting orbit into the satellite orbit. The problem then 
remains of programming the direction of thrust during the first phase to achieve 
maximum fuel economy. Equations are available, from which, by numerical 
integration the solution to this problem may be found"*’. 


3. Programming for Maximum Range of a Rocket Missile 


To illustrate the kind of mathematical techniques which are usefully employed 
when considering the type of problem described at the end of Section 2, a 
method is described of calculating the trajectory of a rocket missile so that it will 
achieve its maximum range over the Earth’s surface. 


Theoretically, optimisation should be carried out with respect to both the 
magnitude of the motor thrust and its direction. However, for a large long-range 
rocket weapon, the optimal magnitude of thrust proves to be impracticably large for 
reasons which have already been mentioned and it will therefore be assumed that 
when the motor is operating, it is always developing its maximum thrust. It remains 
to achieve optimal conditions with respect to the direction of thrust. It is also 
assumed that the period of motion divides into two phases. During the first phase, 
the motor operates at full power and drives the missile towards the vertex of its 
trajectory. The motion during the second phase is one of free fall under gravity 
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Fig. 4. 


alone until the atmosphere is re-entered, when the missile may follow a ballistic 
trajectory under the additional force of air resistance or further extend its range by 
gliding towards its point of impact. 


For simplicity, consider the two-dimensional problem of a missile launched 
from a point on the Equator and moving in the plane of this great circle. Let 
Ox,, Ox, (Fig. 4) be rectangular axes through the Earth’s centre, rotating with it 
and lying in the Equatorial plane and let (x,, x,) be the co-ordinates of the rocket at 
time t. If c is the jet velocity and M is the rocket’s mass at any instant, the accelera- 
tion due to the motor thrust is —cdM/(Mdft). Since the fuel consumption pro- 
gramme is to be assumed fixed, this acceleration will be a known function of 1, 
say f(t). If (/,, /,) are the direction cosines of the direction of thrust during the first 
phase, the equations of motion of the rocket are 


%+Fi=fl, (i=1, 2) 


where (—F,, —F,) are the components of the missile’s acceleration due to the air 
resistance, the gravitational attraction and the centrifugal and Coriolis forces 
associated with the rotating frame of reference. If it is assumed that there is no 
alteration in the settings of any aerodynamic control surfaces and that the orientation 
of the rocket relative to the direction of motion remains invariant during the time 
that air resistance is operative, it is permissible to write 


explicitly. 


Suppose that the functions /,(t) determine the thrust direction programme 
resulting in the achievement of maximum range. Then equations (3) can be solved 
for x,=x; (0), specifying the optimal trajectory. At t=0 the rocket is launched from 
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L ig. 4) and x;, x; take given values. Suppose that the motor operates until t=T, 
“all-burnt” occurring at A. Let A be the point (q,, g.) and let (g,, g,) be the 
velocity components at this point. From A to the point of impact /, the rocket falls 
under gravity. The position of J is completely determined by the-design of the 
missile, the position of A and the missile’s velocity at this point. Thus the range R 
along the Equator from the position of L to the point J is an explicit function of 
the variables q:, 


ie. R=R(q,, Go» Gi» 


This will be a known function of q;, g; for any particular missile. The thrust 
direction programme must now be selected so that q;, g; take values which 
maximise R. 


Consider a small variation in the thrust direction programme in which the 
functions J; are replaced by /;+8/;. Since 1,7+/,2=1, the 4; are related by 
the equation 


but are otherwise arbitrary (the repeated index summation convention is being used 
in equation (6) as elsewhere in the argument). Taking the first variation of the 
equations of motion (3), we obtain 


OF, 


4 oF 


Since the launching point and the missile’s initial velocity are not subject to 
variation, 5x,=8x,;=0 at t=0. Equations (7) accordingly fix the functions 5x; in 
terms of the 5/; and hence the varied trajectory. 


Putting 5x;=5p,, equations (7) are replaced by the first order system 


OF 


8x,—8p,=0. 


Solving these equations by the method of the variation of parameters’, leads to 
consideration of the system of equations 


oF; OF; 
= — 


0 


yi—2z=0, 
for the four functions y,, z; (i=1, 2). Let 


Zit 


G@=1, 2; k=1, 2, 3, 4) 


represent a linearly independent set of solutions of the equations (9). 
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Also, let A be the determinant of the fourth order given by 


A= Z11 Zi2 Zi4 
223 
Mis 


(10) 


Then equations (8) can be integrated over the time interval (0, 7) into the forms 


| 
0 


sp. (1)=84.= | dude, 


where «,;, 8,; are the determinants formed from A by replacement of its i* row by 
the elements 


(T), (T), Yrs (T), Yra (T), 
and Ze Zeg(T), 25 Zee (TF), 


respectively. 


The variation in the thrust direction programme results in a variation in R 
(equation (5) ) of 


oR 
sR= 8q-+ 2G. 8g, 
T 
OR 
0 


after using equations (11). 


If the original value of R was a maximum, 6R=0 for all functions 4/; satisfying 
equation (6). This will be so if 


(aR, 


(13) 


This equation determines the optimum thrust direction programme. 


Integration of equations (3) and (13), determining the optimal trajectory, will 
have to be performed numerically. The presence in equation (13) of quantities 
which have to be evaluated at t=T necessitates performing this integration in a 
sense opposite to that in which the rocket describes its trajectory. The conditions 
existing at “all-burnt” will have to be chosen more or less arbitrarily in the first 
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instance and, after integration back to the launching point L, adjusted until the 
conditions at L are in agreement with those specified. The computing problem is 
clearly one in which it will be necessary to acquire the services of an electronic 
computor. If, however, the problem is simplified by neglecting air resistance, the 
curvature and rotation of the Earth and variations in gravity with height, an 
analytical solution is easy to obtain. Such a solution can be used to assess the 
increase in range to. be expected when any thrust direction programme is modified 
to follow an optimum schedule. This case is considered in the next section. 


4. Maximum Range over an Airless Flat Earth 


In this case it is convenient to take the origin O at the point of projection and 
the axes Ox,, Ox, horizontally and vertically respectively. The missile follows a 
parabolic track from A to the point of impact J where its trajectory intersects Ox,. 
It will be found that equation (5) now takes the form 


g being the constant acceleration due to gravity. 
The equations of motion (3) can be written 
from which it follows that the equations (9) are in this case 


Let y(t), Z(t) be a set of solutions of these equations satisfying the initial 
conditions 


8, being zero unless the two subscripts are identical, when its value is unity. This 
set of solutions is fundamental and from it the determinant A can be constructed in 
the form 


A=}1 0 0 0 
0 1 0 0 
(18) 
0 t 0 1 
It now follows that 
a,,=T-t, %,=0, 
—t, 19 
B= 1, B,,=0, ( ) 
B,,=0, i. 
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Substituting from equations (14) and (19) into equation (13), it is found that 
the optimum thrust direction programme is determined by the equation 


(20) 
+ G2? +2842) 

It will be observed that the thrust direction is independent of the time 1. 


First consider the case when the acceleration f due to the motor thrust 
is constant during the first phase of the motion. Assuming that initially 
x, =X,=%X,=x,=0, equations (15) may be integrated to yield 


x, ft, x,=(L,f-g)t, 
. . . . 
I, and 1, being given by the equations (20). When t=T, x,=q,, X,=q2, X,;=4). 


x,=q., thus providing four equations for the determination of the values of q,, q,. 
Gi» 72. It will be found that, if k is the positive root of the equation 


then qi: =49T? (k+1)3, 
49T*k, 
(k+ 
Ga= gTk. 


Equations (21), determining the trajectory while the motors are operating, may 
now be shown to take the form 


These are parametric equations of a straight: line through O and at an angle @ to 
the horizontal given by the equation 


tand=(1+ 


If ¢ is the angle made by the direction of the thrust with the horizontal, it may 
also be proved that 


4 
tang= 2 = (144) . 
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45° 


Fig. 5. 


It will be observed that as f—> 00, i.e. as the thrust becomes more powerful 
and gunnery conditions are approached, k—> oo and hence @ and ¢ both tend to 
a limiting value of 45°, as might have been expected. 


From equation (14), it is found that 


Complete trajectories for the cases f/g=2, 3, 6, 00, have been drawn in Fig. 5. 
For each trajectory T was selected so that the fuel consumption was the same 
in all cases. 


Equations (25) and (26) show that optimum conditions are achieved by deflect- 
ing the direction of thrust from the direction of motion by an angle 


cot-? [kt (k + 1)? + k#(k+ 1)*] (8) 


in the upwards direction. It is of interest to discover the effectiveness of this device. 
If the thrust magnitude programme is maintained unaltered, but the thrust direction 
is supposed to be always in the direction of motion, the equations of the trajectory 
will be those given in Ref. 15. Assuming that the missile is initially stationary it 
must first ascend vertically. However, there are a single infinity of possible 
trajectories issuing from the launching point and, for comparison with the optimum 
trajectory, we must select that which results in the greatest range being achieved. In 
a particular case it was assumed that f=3g, g=32-2 ft./sec.2, T=50 sec. The 
resulting maximum range, when the thrust was not deflected from the direction of 
motion, proved to be 97 miles. The corresponding solution of equation (22) was 
found to be k=1-271. Equation (27) then yielded Rnsx=104 miles, an increment 
of about 7 per cent. Smaller values of f/g result in a larger percentage improve- 
ment in range, and larger values in a smaller improvement. 


Attention is now directed to the problem of attaining maximum range with a 
rocket whose thrust remains constant. Since the mass will diminish as fuel is 
consumed, it follows that the acceleration f due to the motor will steadily increase. 
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If P denotes the thrust, M, the rocket’s initial mass and y the rate at which propellant 
is consumed, » must be constant if P and the jet velocity do not vary, and hence 
the component of the rocket’s acceleration due to P at time t will be given by the 
equation 


M,-pt 7-t’ 


f= (29) 


where a=P/y,7=M,/p. The equations of motion of the missile are accordingly. 


al, al, 


where, for maximum range, /,, /, are specified by equations (20). 


We now proceed as before and calculate that 


(2—b) p*(p— 

[((2—b) p—1], (31) 
q2= aT (p—1), 


p being the positive root of the equation 


log 
a 
d= of (34) 


It now follows from equation (14) that 


The equations of the trajectory during the period of operation of the motors are 


(9, 
x,=p' (p—b¥X(D (36) 
X,= pX (t)—4 0’, 
log 
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The angle of inclination of the thrust to the horizontal is given by the equation 


tang= (1-2), 


The angle made by the direction of motion with the horizontal at time ¢ is given 
by the equation 


los | 

t 7—T 

log — 


This equation shows that the rocket must be launched at an angle @, to the 
horizontal, where 


In the particular case of the V.2 rocket missile, P= 60,000 Ib. wt., M, =28,500 Ib., 
Ib./sec. Hence a=7,028 ft./sec., == 103-6 sec. Taking the duration of the’ 
thrust to be T=70 sec., it is found that b=1-184, d=3-510. The positive root of 
equation (32) is then p=2-658, resulting in a value for Riuax from equation (35) of 
277 miles. 


Substituting appropriate values for the parameters appearing in equations (38) 
and (40), it is calculated that the angles made by the thrust and the initial direction 
of motion with the horizontal are 


o=49° 12, 0,=23° 23’, ‘ . (41) 


respectively. The portion of the trajectory from the launching point to “ all-burnt” 
has been drawn to scale in Fig. 6. Time, in 10 sec. intervals, has been marked off 
along the curve. 
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The maximum range of this weapon is usually stated to be about 220 miles. 
The difference of 57 miles between this and the calculated figure will be accounted 
for to some extent by the effect of atmospheric resistance which has been neglected 
here. Against this must be set the effects of the Earth’s curvature and the reduction 
in gravity with height. Both these factors will tend to increase the range. It is there- 
fore difficult to assess with precision the improvement in range which would result 
if the thrust direction programme of this rocket were converted to optimum. Kooy 
and Uytenbogaart"® estimate that the effect of air resistance on the range does not 
exceed 6 per cent. It accordingly seems reasonable to expect an improvement of 
about 10 per cent., i.e. of the same order as that found in the previous 
numerical case. 
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The Measurement of High Speed Air Velocity 
and Temperature using Sound Wave 
Photography 


D. E. ELLIOTT, B.Sc.(Tech.) 


(Liverpool University*) 


Summary: A transient sound wave is generated by discharging a spark in an 
air stream, the velocity and temperature of which are required. After a time 
interval of about 25 microseconds a shadowgraph photograph of the sound 
wave is taken; a second photograph is obtained after a further known time 
interval of the same order. The two exposures, both on the same negative, 
show the propagation of the sound wave in the air stream and from this the 
Mach number, the true velocity, the local velocity of sound, and hence the 
temperature, can be calculated. 


Using this method, measurements of Mach numbers of the order of 
0-5 gave values between 98 per cent. and 100 per cent. of those calculated 
from pressure measurements. Typical examples of the sound wave photo- 
graphs are shown. With further development along the lines indicated in 
the paper greater accuracy should be possible. The local speed of sound 
was measured to an estimated accuracy of +1-5 per cent. 


Since only a very short time interval is needed to obtain the photographs, 
the method appears promising for investigating explosions of brief duration, 
or dealing with flows of pulsating character. 


1. Introduction 


In high-speed gas-flow analysis the use of pitot and static tubes together with 
shadowgraph, schlieren and interference photographs is well known. Of these 
methods the shadowgraph and schlieren techniques are generally used to give 
qualitative results, although some quantitative work is done. Unfortunately such 
quantitative work not only involves a good deal of mathematical analysis, but 
seldom gives accurate results. The interferometer is capable of determining the 
change of density of the flow, but extremely sensitive and expensive apparatus, used 
with great care, is needed and this technique is not generally used, except for dealing 
with special problems. By far the greatest proportion of the quantitative work that 
has been done has been based on pitot and static pressure readings. 


From a knowledge of the pitot and static pressures at any point the Mach 
number of the flow can be determined, provided that the ratio of the specific heats y 
is known. A further knowledge of the temperature at the point enables the true 
speed to be calculated. In wind tunnel work, it is usual to assume that isentropic 
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expansion takes place from the known inlet conditions down to the static pressure 
at the point and thus the temperature is calculated; the deviation from this assump- 
tion is negligible when using dry air. Thus for steady wind tunnel conditions pitot- 
static readings will give complete information about the flow. However, when 
dealing with a wider variety of problems the use of pitot and static tubes has the 
following limitations : 


(i) 


(ii) 


(iii) 


(iv) 


(v) 


When used in passages of very restricted area, even the finest search tubes 
cause large deviations of flow. 


When dealing with subsonic flows with Mach numbers above about 0°8, 
special slender pointed static tubes must be used to ensure accurate 
results. Such long points, projecting well upstream of the measuring 
plane, often prove inconvenient in cascade testing. 


Pitot readings in ducts where the gas contains a great deal of moisture 
may give rise to large errors unless suitable corrections are made. 


If it is desired to investigate a rapidly fluctuating flow or a flow of brief 
duration such as an explosion, the inherent damping and time lag which 
occurs with pitot and static tubes and their associated recording apparatus 
make such readings of little value. Moreover if, as in many such cases, 
the temperature is varying, there is an added complication of taking 
temperature measurements at the same point and at the same time. 


Extraordinary conditions may occur where the velocity of a gas is 
required, although the actual composition of the gas is not known, or is 
varying with time. In this case the value y is unknown and therefore 
pitot-static readings would be worthless. 


To investigate such special types of problems several new techniques have been 
suggested or are being developed. A few of these are 


(a) 
(b) 
(c) 
(d) 
(e) 


(f) 
(g) 


(h) 
(i) 
(/) 


The use of hot wire anemometers to record temperature fluctuations. 
The use of special pressure pick-ups to record pressure fluctuations. 
The tracking of radio-active tracers or charged particles. 

Measurement of the rate of discharge from high potential electrodes in the 
stream of fluid. 

The breakdown voltage between two electrodes immersed in a fluid 
(Ref. 6). 

Photographing the passage of the “hot spot” from a spark. 
Measurement of the strength of shock waves by shadowgraph methods 
(Ref. 5). 

Measurement of the distortion, caused by a region of variable density, to 
the shadowgraph image of a network of wires (Ref. 5). 

Measurement of the time taken for an ultrasonic pulse to traverse an air 
stream (Ref. 7). 

The superposition of sound waves on shadowgraph or schlieren 
photographs. 


The present paper deals with the last of these methods. 
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NOTATION 


speed of sound at test section of wind tunnel 


‘a, speed of sound at upstream stagnation conditions 
u___ true air speed 

Mach number 

R _ gas constant 

T temperature 
t time interval 

r, radius of sound wave at a given time 

r, radius of sound wave t¢ seconds later 

d__ distance between the centres of two circular sound waves 
y ratio of specific heats 

8 Mach angle 


2. Proposed Method 

The discharge of a spark in a gas causes a disturbance which propagates as a 
wave of variable density; this wave can be photographed using shadowgraph or 
schlieren techniques. If the disturbance is small the wave will propagate relative 
to the gas with a speed a, equal to that of sound in the region, and after t seconds 
the position of the wave will be as shown diagrammatically in Fig. 1. The sound 
wave will have travelled outwards from its centre C by a distance at feet, and will 
have been displaced downstream from the point of origin O by a distance ut feet. 
Thus, without a knowledge of the time interval, the ratio u/a or the Mach number 
can be found, being equal to d/r where r is the radius of the sound wave and d is 
the distance of its centre from the point of origin. 


4d 


Sound waves in air (subsonic). 


Single exposure. Double exposure. 
M=0°5. M=0°5. 
Fig. 1. Fig. 2. 
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Fig. 3. 
Sound wave in supersonic flow. 
Double exposure. M=2. 


Shadowgraph photographs of sound waves caused by sparks discharging in 
still air were first published by Foley and Souder“). In later work they measured 
the speed of propagation of the sound waves and showed that the initial speed of 
the waves was of the order of twice the speed of sound in the undisturbed air. This 
abnormal speed rapidly decreased as the wave proceeded, and at approximately 
2 cm. from the source the speed was very nearly equal to that of sound. 


By firing a spark in a gas moving with high speed, and taking a single photo- 
graph after a known time interval, the true velocity u of the gas could be calculated 
and would be equal to d/t, where d is the distance between the centre C of the 
sound wave and the origin O, and ¢ is the time interval. However, the condition 
M=d/r would not hold, since the initial speed of the sound wave would be greater 
than the normal speed of sound in the region and the distance travelled by the 
wave in the first t seconds would be more than at feet. For this reason it would 
only be possible to make a rough estimate of the Mach number from such a 
photograph. However, because after a short time interval the speed of the wave 
drops very nearly to that of sound, if two positions of the same sound wave are 
known, as depicted in Fig. 2, the Mach number can be calculated provided that the 
inner ring is outside the region of anomalous propagation. Thus if r, and r, are 
the radii of the wave rings and d is the distance between their centres, 


d 
A knowledge of the time interval t between the two positions enables the true 
speed and the velocity of sound to be calculated, since 


d 
— and a= 
t 
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LIGHT SPARK LIGHT SPARK PIN HOLE (P) SCREEN 
GaP (tz) SHADOWGRAPH 
AIR FLOW 
LENS LENS 
(M2) 
WIND 
TUNNEL SOUND 
SOUND GaAP(S,) 
~ 
t OPTICAL BENCH 
RIGGER ELECTRODE 
OS EARTH ELECTRODE 
HIGH VOLTAGE ELECTRODE 
Fig. 4. 


Projection apparatus (not to scale). 


If the flow is supersonic the sound wave pattern takes the form shown in 
Fig. 3. The previous relationships still hold, but in addition 


Having found a, the temperature at the point can be determined from a°=gyRT 
If the variation of y with temperature is known, any one particular value of a 
completely determines the corresponding value of T. 


3. Experimental Work 


Because of the possibilities that would be opened up by an apparatus working ~ 
on these lines, it was decided to carry out some experiments to see if such a method 
were feasible and if so to determine what accuracy could be obtained. So that the 
apparatus should have wide application it was thought preferable to employ a 
shadowgraph technique rather than a schlieren system, which would need expensive, 
high quality lenses or mirrors and relatively complicated optical benches. 


3.1. MEASUREMENT OF TIME INTERVAL 


To measure the time interval a second spark gap was placed in the still air 
outside the tunnel (see Figs. 4 and 6). This gap was connected in series with the 
wind tunnel spark gap and the light source threw a shadow of the sound waves from 
both gaps on to the screen. By assuming that the wave from the second gap travels 
with the speed of sound in the known atmospheric conditions, the time interval can 
be estimated. It was realised that this method of checking one fact by another of 
the same type is not scientifically rigorous, but the method has the advantages that 
no further apparatus is required, and a single negative can be used to record the 
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progress of both the primary and the timing sound waves. The time interval 
estimated in this way will tend to be slightly longer than the true time interval, 
since the sound wave in the still air will be travelling at a speed very slightly greater 
than the true sonic speed. Thus the true air speed calculated will be slightly low. 


3.2. PROJECTION APPARATUS 


The general layout of the projection apparatus used is shown diagrammatically 
in Fig. 4. The sound wave spark gap S, was situated across a wind tunnel specially 
constructed for these tests, and the timing spark gap S, was placed in line with S, 
but immediately outside the tunnel. 


Two light gaps L, and L, were mounted, independently from the remaining 
electrical apparatus, on a simple optical bench. Lenses M, and M,, of about 6 cm. 
focal length, were securely mounted on the bench so that M, focused the light from 
L, on to the gap L,, while M, collected the light from both gaps and concentrated 
it on a pin hole P (which became a secondary point source for both light gaps). In 
this way a rigid assembly was obtained which could be orientated freely to project 
the beam of light in any required direction. 


The two sound wave gaps S, and S, were fired together (being in series), and 
about 25 microseconds later L, was triggered; L, followed after a further interval 
of about 25 microseconds. A screen or photographic plate N received the light from 
the pin hole and a shadowgraph of the tunnel and the sound waves was formed. The 
distance from the point source to the centre of the tunnel was about 38 in. For 
visual observation a translucent screen was used and it was found that the waves 
could most easily be seen when the screen was about 3 ft. away from the tunnel. For 
photographic recording, the negatives were placed only 8 in. from the tunnel, in 
order to obtain sharp images. 


3.3. ELECTRICAL APPARATUS 


The electrical circuit shown in Fig. 5 can be broken up into three units: — — 

(1) Power Pack. The two metal rectifiers m, and m, were connected in 
opposite directions to one side of the 10 kV. condenser C,, the other terminal of 
which was connected to the 3,500 kV. winding of a 230 V. mains transformer. The 
output from the power pack was supplied through the limiting resistance r, to the 
remainder of the electrical apparatus. Since the power pack was used only to charge 
condensers and no continuous current was drawn, the voltage across the output 
terminals rose to nearly twice the value of the peak voltage of the secondary 
winding. Thus the theoretical voltage=2 x 3,500./2=10 kV. The actual output 
was about 9 kV. 


(2) Impulse Generator. Preliminary experiments showed that the 10 kV. 
voltage, as supplied from the power pack, was insufficient to obtain good 
sparking characteristics across the two sound gaps in series and therefore it was 
decided to boost the voltage by incorporating a two-stage impulse generator. Con- 
densers C,, C, and C, were charged in parallel to 10 kV. through a resistance 
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Fig. 5. 
Power pack impulse generator and delay circuits. 


network r, to r,;, and were discharged in series through the three-ball gap G,, the 
two-ball gaps G, and G,, the resistance r,, and the sound wave gaps S, and S,. This 
discharge of the three condensers in series, giving a total theoretical voltage of 
30 kV., was initiated by applying a pulse from an induction coil circuit to the centre 
ball of.the three-ball gap G,. The resistance r, was used to cut down the illumina- 
tion from the sound gaps S, and S, and so prevent fogging of the photographic 
plates. 


(3) Time Delay. The firing of the three-ball gap G, allowed the condenser C, 
to discharge and thus the voltage at point F dropped to zero. Condenser C,, being 
unable to discharge rapidly, caused a negative pulse to propagate along the - 
inductance-capacity transmission line h,, C;, h., Cy, hs, Cio, A, to the centre ball of 
a light gap. Thus L, was triggered after a time delay, the length of which depended 
on the values of the inductances and capacities in the transmission line. Condenser 
C,, discharged through the light gap L,, providing illumination for the first exposure, 
while condenser C,, caused a further pulse to be sent along the second delay line 
to the light gap L,, which gave the second exposure. 


After preliminary settings of the various ball gaps which consisted of 3 in. 
diameter phosphor bronze balls, consistent results were obtained, the time intervals 
showing little variation. 


3.4. PHOTOGRAPHIC TECHNIQUE 


Many experimenters using spark discharges have found that orthochromatic 
emulsions were just as sensitive to extremely short exposures as were panchromatic 
plates. The author found, however, that he could obtain better results using super- 
speed panchromatic plates. These plates were developed for about 16 minutes in 
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INLET FLARE 


~ 
B 
S44 
| 
TUNNEL SPARK GAP 
STILE AIR SPARK CAP 
S 
' 


— 


Fig. 6. 
Wind tunnel. 


normal developer. This combination proved satisfactory but it is not suggested as 
the best compromise, since very little experimenting was done to determine the 
effects of different developing techniques on different grades of emulsions. It may 
well be that the ball sparks used were much richer in red light than those of previous 
experimenters and hence gave better response with panchromatic materials. How- 
ever, since the red radiations from sparks tend to last longer than the shorter wave- 
lengths, it would be better to use only the blue light, since the sound wave would 
tend to become sharper as the length of exposure decreased. 


3.5. WIND TUNNEL 


Figure 6 is a sketch of the wind tunnel. The side pieces B were machined 
to give a slightly diverging working section, approximately 2 in. wide by 4 in. long 
and 5/16 in. thick, followed by a smooth contraction to a throat 14 in. wide. This 
use of a throat after the working section ensured a constant Mach number in the 
working section, so long as the pressure ratio between inlet and exit was above the 
critical. The divergence along the working section was adjusted to offset the 
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Fig. 7. 
(Magnification = 3). 
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thickening of the boundary layer and hence to give uniform speed over the length 
of the working section. An } in. diameter hole was drilled in each of the glass 
walls and a length of 30 gauge copper wire was sealed into each hole to form a 
spark electrode. Care was taken to leave the inner surface of the walls smooth 
and flush. 


4. Results 


Figure 7 shows a typical photograph obtained when gap S, was short-circuited, 
leaving only S, in the tunnel to produce a sound wave. The dense black disc 
roughly denotes the point at which the sound wave was initiated, and is the shadow 
of the sealing wax which was used to fill the holes through which the electrodes 
were inserted. The dense black lines are the shadows of the conductor wires leading 
to the electrodes. The sound waves are the non-concentric dark and light rings, at 
the centre of each of which can be seen a small irregular dark circle, which is the 
shadow of the hot spot caused by the local heating of the air by the spark. From 
this plate only the Mach number of the flow can be estimated. 


Figure 8 shows a typical photograph obtained when both sound gaps S, and S, 
were fired. This is similar to Fig. 7 with the addition of two concentric dark and 
light rings which are shadows of the sound wave emerging from the still air spark 
gap. These concentric rings enable the time interval between the exposures to be 
estimated if the still air (room) temperature is known. 


In all, about 20 plates were obtained at a tunnel speed of 0:47M as calculated 
from the average of static pressure readings. The plates were analysed by throwing 
an enlarged image on to a sheet of drawing paper and tracing the boundary between 
the dark and light rings of the sound waves. This boundary represents the point on 
the pressure wave curve where the density gradient changes from a positive to a 
negative sign, and it approximates very closely to the peak of the pressure wave. 
In calculating the results account was taken of the different magnifications of the 
wind tunnel and still air sound waves due to their not being in the same plane. No 
other correction is needed with this projection system which gives undistorted 
magnification of the object. 


The values of M obtained from analysis of the photographs varied from 0-46 
to 0-47. Thus the error in measurement as compared with the static pressure 
calculation lay between 0 per cent. and —2 per cent., the negative sign showing that 
all the results were on the low side. 


The value of a/a,, the ratio of the speed of sound in the tunnel to that in the 
upstream stagnation conditions, was found to be within — 1-0 per cent. to + 1-5 per 
cent. of the theoretical value. With an error of 1-5 per cent. in estimating a/a,, 
the error in the value of the temperature calculated from the relationship 
a/a,= /(T/T,) (since y and R are constant), would be approximately 3 per cent. 
(almost 9°C.). The value of the true air speed showed an error of between zero and 
—2°5 per cent. 
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No comparison between the direction of flow of the air, as obtained from the 
photographs, and the axis of the tunnel was made. Inspection of the photographs 
suggests that the direction so obtained would be at least as accurate as that obtained 
by conventional pitot yawmeters. 


During the tests one of the negatives was accidentally given two double sets of 
exposures. The four resulting sound wave shadowgraphs were rather faint. With 
further research it is quite feasible that four or more exposures of the same sound 
wave could be taken without reducing the contrast below an acceptable level. There 
would be obvious advantages in doing this. 


The accuracy obtained in the present experiments does not represent the best 
that could be achieved by this technique. An analysis of the probable causes of 
errors is given in Section 5, together with some methods for their reduction. 


5. Errors 


The errors inherent in this method can be divided into the following 
categories : — 

(1) The discrepancy which occurs between the actual velocity of propagation 
and the true sonic velocity assumed in the calculations. The actual velocity is bound 
to be higher than that of sound in the region by an amount which depends 
upon the strength of the sound spark and the length of time which elapses between 
the initiation of the disturbance and the exposure of the first shadowgraph. The 
energy input to the sound spark should therefore be kept to a minimum, and the 
time interval should be as long as possible. This discrepancy may account for all 
the observed results of Mach number being less than those calculated from 
pressure readings. 


(2) In the theoretical calculations it is assumed that the tunnel flow is uniform 
over the area covered by the sound wave. In practice this may not be so and, in 
order to keep such errors to a minimum, this area, and hence the total time interval, 
should be kept as small as possible. This requirement is in direct conflict with that 
needed to minimise the error dealt with in (1). To overcome this problem and to 
obtain greater accuracy an alternative spark arrangement could be used, as 
follows : — 


Two sound gaps are placed in the air stream separated by a distance of about 
2 in. They are fired simultaneously (although a slight difference in firing times 
would not impair accuracy) and two shadowgraphs are taken as before. Fig. 9 
shows the theoretical position of the resulting sound waves, the value of M being 
given by the ratio (s—/)/(s+J). This method has the threefold advantages that 
(a) The initial time delay can be made much longer, as required by para- 
graph (1). 
(b) The second time interval can also be made longer, and hence the distances 
s and | between the sound waves are increased, giving greater measuring 


accuracy. 
(c) Even with these increases in time intervals, the area over which constant 
flow has to be assumed is reduced. 
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Fig. 9. 
Double spark arrangement. 


A further extension would be to use several gaps at intervals along a tunnel so 
that one double exposure photograph would enable the velocity at several points to 
be established simultaneously. 


(3) Errors due to shadowgraph technique can be reduced by using better point 
sources. With gaps of the Liebessart type projection distances of up to 20 ft., as 
compared with 3 ft. in these experiments, could be used, thereby giving clearer 
sound wave shadows. Alternatively an optical system giving parallel light would 
give sharper pictures. 


(4) Some error is due to the fact that the sound spark line is not always 
exactly straight between the electrodes. In the experiments the electrodes were 
30 gauge (0-012 in.) copper wire cut normally to their length, no attempt being made 
to provide needle-sharp points, since it was thought that needle points would wear 
and that blunt ends would suffice for the experimental period, during which the 
spark gaps would be fired many times. It is estimated that when the photographs 
shown were taken, the sound gaps had already been fired between 200 to 300 times, 
mostly during visual observation of the sound waves. By watching the sound spark 
gaps it could be seen that the spark path was sometimes bent and irregular, this 
being more pronounced with the tunnel gap. This may account for some of the 
woolliness in parts of the sound wave photographs. The use of needle-sharp 
electrodes would greatly intensify the potential field, thus ensuring much straighter 
spark channels, and would be worth while for accurate work. 


(5) The time duration of the sound sparks and the light sources will materially 
affect the sharpness of the sound wave shadowgraphs. No measurements were 
taken of the actual time occupied by these sparks but, from the photographs 
obtained, this duration must have been quite small, probably of the order of 
1 microsecond. Further research could no doubt reduce the duration of both sets 
of sparks and the clearness of the shadowgraphs would be increased. 


(6) The error in measuring the time interval could be reduced appreciably by 
using better electronic time delays, as mentioned later. 
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6. Previous Work 


The double-exposure method is only one of many ways in which sound waves 
could be used to determine air velocities. However, the only other work along 
these lines known to the author is that done by Eisenstein, Clark and Carlson®*. 


In their first series of experiments, sound waves, emitted by a group of five 
sparks placed flush with the walls of a high speed wind tunnel, were photographed 
using schiieren technique. By obtaining several photographs, each with a different 
time interval, and superimposing all the sound waves on one tracing, it was possible, 
by measuring the tangent angles at any point, to calculate the component velocities 
in any desired direction. The errors involved in the angle measurements were large 
and even with a perfect set of tracings, constructed with compasses, the error was 
of the order of +8 per cent. The actual error in the results from the photographic 
tracings was about +15 per cent. 


A second series of experiments was conducted to try to calibrate a section of a 
wind tunnel. To do this, a single spark gap was placed in the nose of a vane placed 
downstream of the region concerned. By obtaining schlieren photographs at 
different time intervals and plotting a graph of distance travelled by the sound wave 
against time, the velocity of the flow and the sonic velocity in the region could be 
estimated within +1 per cent. 


Unfortunately, they were unable to obtain more than one schlieren photograph 
per negative and, since slight fluctuations in flow were inevitable, the different photo- 
graphs were liable to be exposed at slightly different speeds. 


The time delay circuit used by Eisenstein was of a much superior design to 
that used by the author. The variation of the time delay was effected by changing 
the capacity of a capacity-resistance circuit. Theory indicated that the time delay 
would be directly proportional to the capacity and this was found to be completely 
justified when calibrated against oscillographic measurements. Since the capacity 
box could be adjusted in steps of 0-0001 » F. corresponding to time delay increments 
of about 3 microseconds, over a range from 0-0001 »F. to 0-011 uF. a very wide 
and accurate time delay range was available. 


7. Further Developments 


Of many possible ways of improving this photographic technique of flow 
measurement, one of the most promising would appear to be the use of a moving 
image converter. These newly developed instruments can take repetitive exposures 
of down to one microsecond duration and could be used to take several photographs 
of the successive positions of one or more sound waves. In this way a great deal 
more information covering a wider area of flow could be obtained from each 
composite picture. 


*This work is not widely distributed in the United Kingdom. If difficulty is found in obtaining 
this reference, examples of the types of photograph obtained on the same apparatus are 
given in Ref. 1. 
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8. Conclusions 


The experiments have shown that a double-exposure photograph of a sound 
wave proceeding from a spark can be used to measure the Mach number, true speed, 
speed of sound and the temperature of a gas flow. The error involved was of the 
order of 2 per cent. but with further development better results should te possible. 


The method is of great interest since it is fundamental, being independent of 
pressure and temperature measurements. As such it is applicable even in flows of 
non-adiabatic nature, as when a combustible gas is burning during an expansion 
process or in a flow where the composition of the gas is uncertain. 


Since the readings are practically instantaneous, occupying only about 
25 microseconds, the method is particularly promising for investigating pulsating or 
explosion phenomena. 
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Normal Loading on a Wedge-Shaped Plate 


D. E. R. GODFREY 
(Battersea Polytechnic) 


SuMMARY: The equations of thin plate theory are expressed in polar 

co-ordinates and transformed using the Mellin transform. Problems involv- 

ing discontinuous and isolated normal loadings may then be solved in the 
case of the built-in or freely supported wedge-shaped boundary. 


1. Introduction 


A treatment of thin plate theory has been provided by Stevenson"? in which the 
complex variable together with certain combinations of the stress resukants and 
stress couples have been brought to bear on the usual theory. This results in a 
simplification of its expression and development. In a later paper’? Stevenson has 
obtained the stress resultants and stress couples in terms of the transverse displace- 
ment w of the mid-plane of the plate, using curvilinear co-ordinates, and from these 
the required equations in polar co-ordinates may be obtained. The Mellin transform 
has been used by Uflyand® in connection with an isolated normal force on a wedge- 
shaped plate, a paper in which an expression for the deflection is given. 
Woinowsky-Kreiger™ used solutions of the biharmonic equation to solve the same 
problem. The present author has shown? that advantage accrues from the com- 
bination of the Mellin transform, as used originally by Tranter“, and the complex 
variable approach, in the case of generalised plane stress in a thin plate. In this 
paper it is shown that the equations of thin plate theory respond very well to the 
same treatment and at the same time a complete formal solution for the stress 
resultants and couples is presented. 


NOTATION 


It is proposed to follow the specification of stress set up by Stevenson in which 
the stress resultants and couples as obtained by A. E. H. Love in his Theory of 
Elasticity are also averaged over the thickness 2h of the plate. This allows the 


development of a compact notation for the stress couples on similar lines to 


the Pearson notation pg for an average stress component. 


pq=component in the direction of a unit vector q of the average stress resultant 
acting on an element whose normal is in the direction of a unit vector p. 


pq=component in the direction of a unit vector q of the average stress couple 
acting on an element whose normal is in the direction of a unit vector p. 


It will be seen that if pq then pg is a bending couple, whereas pp is a 
torsional couple. 
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% 60 Gr 


Fig. 1. 
In polar co-ordinates these stress couples may best be illustrated on a diagram 
(Fig. 1) in which the plane of the plate is that of the paper. 


Taking the z-axis perpendicular to the plate, there will also be resultant shear 
components Fz and 


It will be found convenient to put the stress couple components in complex 
combinations in the form 
I¥=r0+6r+2i00, 


where i= /(—1) and (r, 4) are polar co-ordinates. These are related to the corre- 
sponding Cartesian combinations 


T= xy + 9x +2iyy, A= xy yx « « 
We shall also write W’ = rz +i0z=e- (xz + iz). 


Other notation is 


r,@ polar co-ordinates in the mid-plane of the plate 
z Cartesian co-ordinate perpendicular to the mid-plane 
Poisson’s ratio 
E  Young’s modulus 
q normal pressure on the face z=h of the plate 


Note—A bar over a symbol denotes the conjugate complex quantity. A double 
bar denotes the Mellin transform with respect to r. 


2. The Mellin Transform and its Application to Plate Theory 


2.1. THIN PLATE THEORY 


If a plate is loaded normally with a pressure q(r,4) the deflection w of the 
mid-plane satisfies the equation 


Vw= 9) (5) 
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where D, the flexural rigidity of the plate, is given by 


2Eh*® 


Also P=D/(2h) is used as a plate modulus. 


The components of the stress resultants and couples may then be expressed 
in terms of w by the equations 


ro=(1— yp {t (8) 


Forming now the expressions mentioned in equations (1) and (4), 


A’= —(1+v)PV?w : ‘ 
n= 
[ —2i { r (Aw) — Aw \] (13) 
i 


0 
where A 56° 


2.2. THE MELLIN TRANSFORM 


The Mellin Transform f(p) of the function f(r) is defined by 
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Transforms of derivatives of f(r) may be obtained, on integrating by parts, in 
the form 


fr (r) dr=(-— 1% p(p+))...... (p+n— 1)F@). n=1,2,3,4. (16) 


The equations of thin plate theory may now be transformed when it is found from 
equation (5) that 


and from equations (12), (13) and (14) that 
(PA)= —(1+») P(A?+ p?)w 
(FT) . 9) 
(PW) = —iP(A?+p?){A+i(pt+2}w. . .  . (20) 


The success of the method is due to the fact that the differential equation (17) 
is easily capable of solution and its complementary function will be written in 
the form 


Ae + Ae- 64 Bel +264 


since w will be real if p is real. The particular integral will assume various forms 
according to the nature of the loading. 


2.3. BOUNDARY CONDITIONS 
Consider an infinite wedge-shaped plate bounded by 6= +2, 


(a) with built-in edges, giving w=0 and Aw=0 so that the transform satisfies 
w=0, Aw=0, . . . . . (22) 


(b) with freely supported edges, giving w=0 and @r=0 which, on using the 
fact that 0w/dr and 0?w/dr? are both zero on the edges, leads to 
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2.4. INVERSION INTEGRAL 


Provided that | r*-1 | f(r)| dr is bounded for some k > 0, an inversion integral 
0 


may be expressed in the form 
1 
f= dp 


where c>k>0. 


(24) 


Once w has been determined the stress resultants and couples may then be 


obtained in the form 


c+ico 
(1+yv)P 


A’= 


c—too 
c+ico 


c+ico 
Ff dp. 


e—tao 


3. Uniform Loading over a Sectorial Area 


For a constant loading gq we have 
where p>-—4. 
Using equation (17), w will take the form 


w=Ae?+ 4 Bel + 


D 


so that for built-in boundaries 6= +2, condition (22) gives 


qd?** sin (p+2)a 


A=~ DO+4) 

sin pa 

2DH (p,2) p(p+2)? (p+ 4) 
where H (p, 2)=sin (2p +2)2+(p+1)sin2z2. 


(25) 


(26) 


(27) 


(28) 


(29) 
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From these constants the stress resultants and couples may be found at any 
point of the wedge, using relations (25). In particular it is here proposed to find 
the fixing couple 6r in the case of the semi-infinite plate, i.e. when 2=2/2. This 


will necessitate finding both A’ and I”. 


Now (A + p*)w =8B(p+ + when 


qd? +4 1 
p(p+2)(p+4)’ 


so that, from (25), 


(14 


(30) 


Again 


qd’+* 1 
p(p+2)(p+4) 


(A+ip){A+i(p+2)} w =4(p+ !? [B(p+2)— Ap] — 


1 
p(p+2)(p+4)’ 


c+iao 


Equations (30) and (31) show that on the built-in boundary 


(1 —v) A’= -(1+)I” 


or 66=0 and ro=-vor. . . G2) 


Equation (32) is in agreement with the more general result obtained by Stevenson, 
in which he shows that if mn and s are unit vectors normal and tangential to a built-in 
boundary, then — vag. 


It follows then that 


e+tao 


This integral may be evaluated directly by taking c >0 and using a semi-circular 
contour of large radius R with centre at p=c, and which encloses the poles p=0, 
—2, —4, giving finally 


~ 16hr? 
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4. Isolated Normal Force at a Point on the Axis of Symmetry Distant 
d from the Vertex 
Ufiyand™ has provided a compact way of representing an isolated normal force 
by considering the Laplace transform of equation (17), which we shall rewrite for 
this purpose as 
(A? + p*){A? + (p+ 2)"} w=f 


For the moment only the particular integral, which will be written ®, (4), is 
considered. 


The Laplace transform of a function F (@) is given by 
co 


so that, utilising the usual machinery of the Laplace transform, we have 


+ (0+ (8)= | ef, (0) 


1 


The convolution theorem then leads to 


6 


A form of ®, (@) corresponding to an isolated force may now be obtained by con- 
sidering the pressure g to act only on a small element of the surface enclosed by 
+8, r=d, r=d+a. 


From equations (17) and (35) 


ad?+* 
=0, e<|0|<a, 
Therefore 
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Letting «—> 0 and a—>0, while keeping gqdas=4W, we obtain 


&, (6)= Wad?*? sin p@ 


8D(p+1) 
_ Wart? fsinp@ sin(p+2)6 
8D(p+1) p pt+2 


(38) 


The solution of equation (35) may now be written 


W=2A cos . . . (39) 


and the constants are determined from the boundary condition, which will be taken 
to correspond with built-in edges. We then have 


p—J(p,2) 


2A= 
Wa*? p+2—K(p,2) 
where J (p, (p+ 1) cos 2z—cos(2p+2)2 
K (p, 2)=(p+ 1) cos 22+ cos (2 p+ 2) z. 


The stress resultants and couples may now be expressed formally by using relations 
(25). Integration for a general value of « will inevitably involve numerical com- 
putation, but for the semi-infinite plate the integrals are evaluated and agreement 
is obtained with known results. In this case 


| 
16D p p+l 


tan p ~ cot p = 
+2 9 
Wd? 2 2 


16D p+2 p+l 


As in Section 3, 6r= — ere — 
~ d dp 


sin p 
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Taking c=0, p=in, and allowing for the pole at p=0, we have 


Wd? sin mn 
2 
0 
where m=log, (d/r). 
Using the result 
. 
sin mn 
dn=tanhm 
sinh 


gives finally the fixing moment 


Wd? 
+d)’ 


which is in agreement with the value given by Woinowsky-Kreiger. 


5. Conclusions 


The Mellin transform seems to provide a natural machinery for solving plate 


problems with a wedge-shaped boundary. Indeed, methods using the solutions of 
the biharmonic equation in polar co-ordinates usually have to be knitted together 
into the form of integrals such as arise from the inversion integrals. The use of 
stress resultants and couples in the form of the combinations A’, I’, ¥’ simplifies 
the formal expression of results. 
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Compressible Subsonic Flow in 
Two-dimensional Channels* 


Part 1: Basic Mathematical Theory t 


L. C. WOODS, A.F.R.Ae.S. 
(Department of Applied Mathematics, University of Sydney) 


SUMMARY: Equations for the calculation of the subsonic flow of an inviscid 
fluid through given two-dimensional channels (the “ direct’ problem), and 
for the design (the “ indirect” problem) of such channels are derived. The 
method is based on conformal mapping, and in the special case of channels 
with walls made from a number of straight sections, or with wall pressure 
prescribed as step-functions, yields the same results as the well-known 
Schwarz-Christoffel mapping theorem technique. However, it is more 
general than this latter method, since it is capable of dealing with curved 
walls or continuously varying wall pressures. The compressibility of the 
fluid is allowed for only approximately, the ideal gas being replaced by 
a Karman-Tsien tangent gas. 


In Part II the theory is applied to various problems of aeronautical 
interest, perhaps the most important of which is to the setting of “ stream- 
lined ” walls about a symmetrical aerofoil placed in the centre of the channel. 


1. Introduction 


It is well known that the flow of an incompressible fluid through channels with 
walls made from a number of straight sections, or with walls, the pressures on 
which are prescribed as step-functions, can be calculated quite simply by making 
use of the Schwarz-Christoffel theorem on the mapping of the interior of a polygon 
on the upper half-plane’. This paper presents a method of dealing with the 
compressible flow between curved walls (the direct problem), or between walls on 
which continuously varying pressures are prescribed (the indirect problem). The 
direct problem is shown to involve the solution of an integral equation, which in 
general requires an iterative method of solution, whereas the indirect problem can 
be solved immediately without iteration. 


The indirect problem has been considered by a number of authors. Lighthill"’, 
in his paper on aerofoil design, also dealt with the design of symmetrical channels, 


*The work described in this paper was done early in 1952 while the author was a memter of 
the New Zealand Scientific Defence Corps, seconded to the Aerodynamics Division of the 
National Physical Laboratory, Teddington. 

+Part II will be published in Novemter 1955. 
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while Stanitz has designed channels, both by a relaxation method™ and by a method 
very similar to that described in this paper“*. 


There are many other papers, too numerous to cite here, dealing with various 
aspects of channel flow, but it seems that few papers have been written that give 
an exact treatment of the direct problem for flow between curved walls, and the 
method given in this paper is probably new. 


In some channel flow problems the boundary conditions are such that the 
pressure is prescribed on some sections of the walls, and the shape is prescribed for 
the remaining sections. These mixed boundary conditions have been considered in 
a recent paper by the present author’, but one special and important case of the 
“mixed problem” is treated in this paper. This is the flow between a constant 
pressure well and a wall whose shape is prescribed. This work has an important 
application to the theory of flexible-wall wind tunnels, and this will be given 
in Part II. 


Compressibility is allowed for approximately by making use of a “ Karman- 
Tsien tangent gas.” This is a gas whose adiabatic curve in the p, 1/p diagram is a 
straight line tangential to the curve for an ideal gas. The point of tangency is taken 
to be the point corresponding to undisturbed stream conditions. The theory is 
exact for such a gas, which includes the ideal incompressible gas as a special case. 


The theory of the tangent gas, appropriately modified for the particular 
application to be made in this paper, is given in the next section. 


NOTATION 


The most important symbols used in this paper are defined by: — 
(x, y) _ the physical or z-plane, where z=x+iy 
n, s distances measured normal to and along a streamline respectively - 
(q. 9) velocity vector in polar co-ordinates | 
local and stagnation densities respectively 


H, K channel widths in the z-plane at infinity upstream and downstream 
respectively 


(9, ¥) plane of velocity equipotentials (¢=constant), and streamlines 
(Y¥=constant) 


U, V _ velocities at infinity upstream and downstream respectively 
Mach number 
B = (1-M?) 
Q — defined by equation (4) 


*The author’s theory for the indirect problem was developed independently of Stanitz’s work, 
and was given in A.R.C, 14,838, April 1952. 
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m_ defined by equation (5) 
f = 
C, _ the pressure coefficient 
9, values of 6 on and Y=h/m,, respectively 
rt the deflection angle of the channel 
h_ the channel “width” in the w-plane, where w=¢ + im,¥ 


Suffix co denotes values in the undisturbed stream at infinity. 


2. The Theory of the Tangent Gas 


We start from the equations’ 


aq _ 
én 
os qon 


' where (q, 9) is the velocity vector in polar co-ordinates, n, s are distances measured 


normal to and along a streamline respectively, and M is the local Mach number. 
They can be transformed into 


20 an 
oy p a) 
p an 
= =0, 
by writing d¢=qas, dy= 
0 
— M2) 


q 


where p, p, are the local and stagnation densities respectively. From (2) it is seen 
that 9, ¥ are the potentiaal and stream functions respectively. Equation (3) yields 


B 


where U is some standard or reference velocity, 8=(1—M7)!, and the constant of 
integration has been selected so that 2=0 when qg=U. 
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equations (1) can be written 


_ 
oy 6 
mop 
1 1/(y-1) 
Now m= (14+ 
where y is the ratio of the specific heats. Hence for M< 1, 
Let the subscript 00 denote values in the undisturbed stream at infinity; then 
Hence the approximation 


due to von Karman’, is reasonable in subsonic flow. 


At this stage we change our point of view and regard equation (7) as an exact 
relation for some imaginary gas. The right hand side of the equation is taken to 
be the value of m,, for an ideal gas, while the left hand side is the value of m 
for some imaginary gas. In the following theory all quantities distinguished by a 
subscript 00 apply to the ideal gas, while the other quantities apply to the imaginary 
gas. 


Equation (7) can be written 


of (fe). « 


a’ 


where a is the speed of sound. Using Bernoulli’s equation, namely, 


dp 
—+qdq=0, 
qaq 
dp 
and 
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in equation (8) it can be shown that 


dp 
2 2 2 


Equation (9) shows that 


d(1/p) 


i.e. the imaginary gas has its adiabatic parallel to the tangent of the ideal gas 
adiabatic at (p., 1/p,.). Integrating (9) gives 


the constant of integration being selected to make the two adiabatics touch at 
(Px» 1/p.). The gas defined by equation (11) is a special case of a “ Karman-Tsien 
tangent gas.” 


A more important property of this gas, from our point of view, follows from 
equations (6) and (7). It is that equations (6) can be written in the form of 
Cauchy-Riemann equations : — 


a. 
O(m.¥) 
showing that the complex function 
f= OQ + id, ‘ (12) 
is an analytic function of 


In incompressible flow Q=log(U/q), m.=1, and then 


f=log (2 =log (22), 


where z=x+iy is the physical plane. The only significant difference between the 
theories for compressible and incompressible flow is that the relation 


1 


does not hold in compressible flow. 
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To calculate a relation corresponding to (15) for compressible flow, p and q 
must be expressed in terms of 2. From equations (8) and (10) 


Thus afi=«,,f.., 
q g_ 
M B= M,, Beis 
q 
and hence i (16) 
Solving (16) for 8 and substituting in (4) it is found that 
= sinh cosech (2 + ¢), (17) 
and B=tanh (CQ + ¢), ‘ (18) 
From (8) and (18) 
p= tanh (+2). 


Now it is easily verified that dz=e"ds, dz=ie dn, so that equations (2) yield 
oie ; Po 
This can be reduced after some algebra to 


a= 


clan}, 


where A is the conjugate complex of A, by making use of equations (5), (17), (19) 
and (20). Equation (21) is the required generalisation of (15), and it shows that 
z is not an analytic function of w except for incompressible flow. 


From equations (11), (18) and (20) it follows that the pressure coefficient, 
C,=(p—px)/ GexU*), 
2 
= {1 — B,, coth (2+ 


i.e. C,=2 sinh 2 cosh cosech (2 + ¢). 
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3. A Theorem on Analytic Functions in an Infinite Strip 

The independent variables of the following theory are 9 and y as defined by (2) 
and not the usual variables x and y of the physical plane. As shown in Section 4, 
this leads to some difficulty in determining the actual boundary values as functions 
of ¢, but on the other hand the shape of the boundaries is considerably simplified, 
for now the infinite strip —co<¢=< 0, OX m,v<A, in the w-plane represents 
flow in a channel of quite general shape. 

The dependent variables of the theory are 2 and 6, where 2 is related to the 
velocity by (17) and (19), or to the pressure coefficient by (22). From (12) and (13) 
it is seen that the problem is now to determine f as a function of w in the infinite 
strip just defined, subject to the given boundary conditions. These boundary con- 
ditions are that for the direct problem the imaginary part of f (@=tan~' (dy/dx)) 
on w=0, w=ih, is known or can be deduced, while for the indirect problem this is 
true of the real part of f(©). A theorem appropriate to this problem will now be 
stated and proved. 

Theorem. If y is an analytic function within the infinite strip —°o< x= 00, 
0<y<A, is finite at infinity, and its boundary values satisfy the 
condition given below, then its value at any point P(z) within the 
strip is given by 


co 


(23) - 
where a and #8 are the real and imaginary parts of y, 
B,=lim B(x, y), B.=lim B(x, y), and 2, 2, are the values of 2 
at x= 00, — 0, The boundary values 8,, 8, must be continuous 
functions, except possibly at a finite number of points where they 
may have either logarithmic singularities or simple discontinuities. 


If such a singularity is at z=x,, for example, then 


y ~ K log (x, —z)=K log | x,—z| +iK arg (x, —2), 


and so if K is real (imaginary), 8, («,) has a simple discontinuity of Kx at z=x,, 
while 2, (8,) has a logarithmic infinity at z=x,. 


Proof. As y is an analytic function within and on the rectangle 
-~-R<x<R,0<y<A, provided it is suitably indented to exclude the logarithmic 
singularities of y, it follows from Cauchy’s theorem when n= +1, +2,..°. , and 
from Cauchy’s integral when n=0 that 


R 
nX+0,0 x+2nih-z x+2nih+ih—-z 
h 
y(R, y) y(—R,y) } 
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Here the radii of the indentations excluding the singularities have been made to 
tend to zero, it being easily verified that the contributions to the integrals from these 
indentations also tend to zero. By addition, equations (24) yield 


R 
X+2nih-—z x+2nih+ih-z “+ 
h 
y(R,y) Ry) | 
n=—N R+2nih+iy—z —R+2nih+iy—z dy. (25) 
0 
lim 3 (5) = (52). (26) 


and in this limit the series in (25) are uniformly convergent as P(z) is not on either 


z=0, or z=ih. Hence lim | = { tim, and, from (25) and (26), 


R 
1 
y(z)= { ¥(x.+0)coth (x, h—O)coth (x-+ih—2) | dx + 
-R 


h 
1 
+ {7(R.y)coth (R +iy—2)-y(-R, yeoth = (-R+iy-a} dy. 


0 


Next we take the limit R —> ©, to find 


y(z)= { 0)coth 2)—y(x,h—O)tanh (2) } AX +4 (Ya + 
(27) 


If the same argument is applied to the conjugate point P(Z), instead of P(z), 
an equation similar to (27) is obtained, but with z replaced by Zz, and, since the 
points P(Z+2inh), n=0, +1, +2,..., all lie outside the rectangle -R<x<R, 
0=<y<A, with the left hand side of the equation replaced by zero. Hence 


co 


0= { ¥(x,+0)coth (x, h—O)tanh (x-2) } AX +4 (Yo 


(28) 


The theorem now follows by adding the complex conjugate of (28) to (27), and 
putting 


y (x, +0)—¥ (x, +0)=2i8,, y (x, h—-0)—y(x, h-0)=2i8,, 


Yoo + 20 = + Y = 2 
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An important relation between the boundary values can be derived from (23) by 
taking either of the limits z—> + ©. It is found that 


| . . . . 


Suppose z has the value x, +i¢ (or x,+ih—ie), then equation (23) still applies 
in the limit « —> 0, but now the Cauchy principal values of the resulting improper 
integrals must be taken. It is not sufficient merely to put z= x, in (23), although 
this does give the correct result for the real part, but the limit z—> x, must be 
properly evaluated. 


4. Compressible Flow Through Channels 
4.1. THE BOUNDARY CONDITIONS 


First it is necessary to assign boundary conditions an infinite distance upstream 
(¢= and downstream (¢=°0). Suppose 


lim (q,4)=(U,0), 
and lim (q,4)=(V, 7), 
then from (4) and (12), 


lim f=0, lim f= log 7) tir. (30) 


¢> 


The value of h, the channel width in the w-plane, can be determined as follows. 
Let H be the channel width in the z-plane upstream at infinity; then from the 
imaginary part of (21) 


14+ 8, 
- 


i.e. h=8,.UH. ‘ ‘ ‘ P ‘ GL 


iH= 


} 


The channel width downstream at infinity, say K, is given by 


from (21) and (30). Hence, using (31), we find 
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4.2. THE Direct PROBLEM 


Equations (23) and (28) can now be applied directly to the function f=0+i0 
in the w-plane (w=9+im,,¥) to give 


f(w)= { 6, coth (¢—w)—6, tanh (9 —w) } dg+4log(Z), (33) 


—co 


the last term of which follows from (30), and 


log (5) =— (0, —9,) d9. (34) 


These are the basic equations in the direct problem, and in the general case their 
solution requires an iterative process, as will be shown. 


First it is convenient to integrate (33) and (34) by parts to obtain 


co co 


=f log sinh (9 — w) (9) + log cosh (9—w)d6, (9)+410g (2), 
(35) 
and )= (9)- | oda @, . .  . (36) 


the integrals in which are Riemann-Stieltjes integrals, to allow for possible dis- 
continuities in #, and 6,. The integrated terms of (35) and (36) vanish since, 
from (34), 

lim (0,—9,)¢=0. 


Now when @, is continuous 


dé, ds, 


ds, do G7) 


d0,= 


where R, (= —ds,/d9,) is the radius of curvature of the wall; a similar equation 
holds for #,. On ¥=0 at w=q’, (35) and (37) yield 


Tu log sinh (¢ +25 Tn log cosh +4 log ), (38) 
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where 79;, are the possible discontinuities in 4,, 6, occurring at ¢=9,, i=1,2,...7n. 
Equations (17) and (38) can now be combined to form an integral equation for gp, 
and a similar one can be deduced for gh by putting w=¢’+ih in (35). Thus we 
have to solve a pair of simultaneous singular integral equations. They are of such 
complexity that in the general case an iterative numerical method of solution 
must be adopted. The procedure is as follows. 


At some point well upstream, where the channel walls are parallel and the 
flow is almost uniform, we draw a normal to the channel walls and take this to be 
the equipotential g= —a, say. Then approximate values of q,(¢) and q, (9) are 
assumed and substituted in 


s= | 


to yield s, s.(¢) approximately. Now R,(s), Ri (s) can be deduced from the 
given wall shape, so that we now have R,, R,, and finally R,qo, Rign as approxi- 
mate functions of ¢. These are substituted in equation (38) and the corresponding 
one for ©, (¢’); from (17) this yields new values of q, (¢), qx (), and the first iteration 
is complete. The ratio of the functions g,(¢), q:(#) initially assumed must be 
adjusted so that the left hand side of (36) agrees with the value derived from (32), 
and this adjustment is necessary after each iteration. The calculation is repeated 
until q,, gn are unchanged by an iteration. It is not proved in this paper that this 
process converges, but “ practical” convergence is shown to occur in an extreme 
case calculated in Part II of the paper. 


4.3. THE INDIRECT PROBLEM 


The equations for channel design are found by applying equations (23) and 
(28) to the function if = —6+i in the w-plane. Corresponding to (33) and (34) we 
find 


{ 0, (—w)—, tanh = (ow) } dothir, . (40) 
co 
On Y=0, w=¢’, equation (40) yields 


and a similar equation holds for 4, (¢’). 
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Now in the design problem we start from given functions q,(s), gn (s) (or 
pressure distributions), and deduce 2, (s), 9,(s) from (17) and (19). Then 9(s) 
for each wall follows from 


| ads, 
(43) 


o= { ands, 


enabling {2, (¢), ,(¢) to be deduced. These values are substituted in (42) and the 
corresponding equation for 9, (¢’), to yield 6,(9’) and 6,(9’). Finally the wall 
co-ordinates are obtained from 


(44) 
_ sind, \ 
Jo 


and similar equations for x,, yn. 


The deflection angle is given by equation (41). If it is desired to specify this 
angle initially, the velocity distributions can be defined only to within an arbitrary 
parameter initially. This parameter is then selected so that (41) is satisfied. 


4.4. THe “CHANNEL” OF INFINITE WIDTH 


Flow past a symmetrical aerofoil at zero incidence placed midway between 
symmetrical walls can be regarded as a channel flow of the type considered in this 
paper if only half the flow pattern is considered. If the channel width tends to 
infinity we have the case of a symmetrical aerofoil at zero incidence in an infinite 
stream. The equations for this case are well known, and are readily derived from 
the more general equations already given by taking the limit h—>0. In this limit 
equations (32), (33), (34), (40) and (41) yield 


6,d if Q,d 
o¢-w «= )¢o-w 
-00 
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5. Further Equations for Incompressible Flow 


The equations given in Section 4 are, of course, directly applicable to incom- 
pressible flow; it is necessary simply to put M,,=0 and 8,,=1 in these equations. 
In addition to these equations others of a similar type can be found for incom- 
pressible flow. These alternative equations are not as useful as those already 
described, but they are worth noting. 


The theorem established in Section 3 is applicable to any analytic function 
which is finite at infinity. Consider, for example, the “displacement function,” 
Uz—w. This analytic function remains finite upstream and downstream at infinity 
if the channel is the same width upstream at infinity as it is downstream at infinity, 
and if the deflection angle + is zero. Suppose this is the case. On ¥=0 and 
y=h/m,, the imaginary parts of Uz—w are Uy, and Uy,—h, and so it can be 
deduced from (23) that 


co 


Uz=w+ 7) tanh dg, 44) 


where the origin of the z-plane is selected to make 


lim (Ux—¢)+ lim (Ux—9)=0. . ‘ . (45) 


From (29) and (45) it is found that y, and y, must satisfy 


2h 


In the limit h—> 00, equations (44) and (46) yield 


o—w 
and lim Ux-9= | vedo. . (48) 


Differentiating (44) and (47) with respect to w gives an equation for 


dz Ucos6 .Usin@ 
= +i 


dw @q q 


which is sometimes useful. 
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Equations similar to these have been used by Thom"? in his “ influence factor” 
method of dealing with channel flows. Goldstein and Lighthill®’ have also used 
similar equations in connection with the Janzen-Rayleigh method of computing 
compressible flows. 


Another analytic function of some value in incompressible flow is 
dw/dz=u-—iv, where (u, v) is the velocity vector in Cartesian co-ordinates. This 
function remains finite at infinity upstream and downstream, provided only that the 
ratio of the channel widths up and down stream is finite, and therefore it can be 
used for channels of more general shape than can the displacement function. 


For this function it is deduced from (23) and (29) that 


u—iv= v, coth —w)—v, tanh (¢--w) ¢ (49) 
co 
1 | 

where Ux — = do. 

In the limit h—> 00 these yield | 

| Us do=0 | 


“ Linearised”” approximations to all these equations can be obtained simply 
by replacing ¢,w by Ux,Uz under the integral signs. This linearisation is valid to 
first order in the magnitude of the boundary values @,, 4, yo, yx—A/U, v, and vp. 


6. Channels in which One Wall is a Constant Pressure Wall 


In an important class of “channel” flow problems, one wall is a constant 
pressure streamline, while the other either has a known shape (most common case) 
or a known pressure distribution. 


The solution for the known pressure distribution follows directly from equation 
(40), since on the constant pressure wall Y=h/m,, say, q,=U, and from (17) 
=0; thus 


co 
f=hir-> 2, coth = (o-w) de, 
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where, from (41), ©, must satisfy 


When @,, and not {2,, is known, the equation corresponding to (51) can be 
obtained as follows. Assume that #,=4, in (33), then there results 


since, from (34), V=U. On w=9’+4ih, this equation becomes 


f= sech 9’) a¢, 


from which it is seen that 2 vanishes on Y=4h/m,,, which is thus a constant 
pressure streamline. Thus if A is replaced by 2h in (53) the resulting equation, 
namely 


fo, cosec (¢—w)d¢, P . (54) 


2h 


clearly applies to flow between a wall of given shape (UV=0), and a constant pressure 


wall (¥=h/m,.). 


Equation (54) yields 
sech (0—9’)d 55 
on ¥=h/m,,, and 
1 co 


ony=0. The first of these gives the shape of the free streamline, while the second, 
in conjunction with (22), yields the pressure distribution on the given wall. 


An alternative derivation of these equations is to be found in Ref. 5, where 
the more general problem of a varying pressure distribution on the “free” stream- 
line is treated. 
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The Effect of Sting Supports on the Base 
Pressure of a Blunt-Based Body in 
a Supersonic Stream 


I. S. DONALDSON 


(Fluid Motion Laboratory, University of Manchester) 


SumMMARY: Experiments have been made to find the effect of the ratio of 
sting to base diameter on the base pressure of an axially symmetric body at 
zero incidence in a supersonic stream. The Mach number of the flow was 
1-994 and the model boundary layer was turbulent. The model used was a 
one inch diameter circular cylinder without boat-tailing. It passed through 
and was supported upstream of the nozzle throat. This method of support 
allowed measurements to be made in the important (and hitherto unexplored) 
case of zero sting diameter. 


As the sting to base diameter ratio was increased from 0 to 0-85, the 
base pressure decreased. The minimum value reached was approximately 
0-8 of the value it would have at the base of a two-dimensional body with a 
similar ratio of boundary layer thickness to base height. The base pressure 
coefficient rose rapidly to zero as the ratio was further increased to unity. 


Under the conditions of the experiments, with a sting to base diameter 
ratio of 0-4 the base pressure coefficient differed from that without a sting by 
approximately ten per cent. With the more modest ratio of 0-2, the difference 

was approximately three per cent. 


NOTATION 


P» base pressure 
p, free stream static pressure 
p, _ free stream density 


u, _ free stream velocity 


C, = ori, . base pressure coefficient 
D _ body diameter 

sting diameter 

h two-dimensional base height 

8 boundary layer thickness 


, boundary layer displacement thickness 


Received February 1955. 
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BASE STING 


MIXING REGION 
Fig. 1. 
Sketch of the flow at the base of an axially symmetric body with a sting support. 


1. Introduction 


Figure 1 shows the flow over an axially symmetric body with a flat base and 
a sting extension. The boundary layer separates at the corner of the base and forms 
an annular mixing region which becomes attached to the sting some distance down- 
stream. As the mixing region approaches the sting it thickens due to mixing and 
due to a decrease in mean radius. The mixing region, the sting, and the base 
enclose a “dead air” region in which the velocities are small compared with the 
free stream velocity, and the pressure is approximately equal to the static pressure 
outside the mixing region. 


The mixing region entrains air from the “dead air” region. However, part of 
the comparatively low velocity air in the mixing region has insufficient momentum 
to overcome the pressure rise through the recompression shock wave and is driven 
back into the “dead air” region. The reverse flow is shown qualitatively by arrows 
in Fig. 1. In the steady state, the Mach number downstream of the expansion fan 
is such that the rate of loss of mass from the “dead air” region by entrainment is 
just balanced by the rate of increase in mass due to the reverse flow. 


When sting supports are fitted, two parameters are important; these are 
(D—d)/D and 6/(D-—d). For small values of d/D, (D—d)/D is nearly equal to 
unity and there is little difference from the flow without a sting. As d/D is 
increased, (D—d)/D decreases but remains, at first, large compared with /(D—d). 
In this range the flow tends to the two-dimensional type. The base pressure tends 
to the two-dimensional value which, for a turbulent boundary layer and small 6/h 
at a Mach number of 2, is almost independent of 5/h. However, the two- 
dimensional value corresponding to (D—d)/D tending to zero is never reached, 
because, as d/D is increased, (D—d)/D and 6/(D-—d) become of the same order of 
magnitude. The length scale perpendicular to the stream direction is then of the 
same order as the length scale in the direction of the stream and the pressure rises 
rapidly to the free stream static pressure. 


With sting-mounted models, strength considerations often limit the ratio of 
sting to base diameter to values greater than about 0-3"). Previously this has 
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Fig. 2. 
Installation of the model in the tunnel. 


precluded the investigation of the flow at the base of a blunt-based body as the 
ratio of sting to base diameter tends to zero. 


Since the experimental arrangement described here does not depend on the 
sting for support, no such limitation is placed upon the sting size, and the flow has 
been investigated for the whole range of values of the ratio d/D, from 0 to 1. 


2. Experiments 


The experiments were made in a vacuum-operated intermittent supersonic wind 
tunnel with a working section 4 in.x5 in. The Mach number was 1-994 and the 
stagnation pressure was approximately atmospheric. 


The flow was observed with a conventional two-mirror schlieren system, and 
photographs were taken with a-spark discharge which had a duration of less than 
One microsecond. 


Base and static pressures were measured with absolute mercury manometers. 
Pitot pressures were measured with a mercury U-tube manometer. 


3. Description of Model 


The model, shown in Fig. 2, consisted of a one inch outside diameter tube 
which passed through the nozzle throat and was supported at its upstream end by a 
vertical tube. The model was also supported vertically at B and horizontally at 
Band C. This method of support made it possible to set the model accurately on 
the tunnel centre line. A cylinder, of 1 in. diameter and 2 in. long, was screwed into 
the horizontal tube at A. 
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Fig. 3. 


Method of supporting the stings. 


Sting extensions were screwed into the centre of the base of the model (Fig. 3). 
The downstream end of each sting fitted into a cross-support in the tunnel diffuser. 
Twelve different stings with diameters ranging from 0 in. to 0-969 in. were used. 
Two 0-015 in. diameter holes were drilled at opposite ends of a diameter of each 
sting, as close as possible to the base of the model. The holes were connected to 
an absolute mercury manometer through the centre of the sting and the base of the 
model. 


4. Effect of Model on the Flow 


Before the main experiments were made, the effect of the body on the Mach 
number and uniformity of the flow in the tunnel was investigated. No stings were 
fitted in these experiments. 


To verify that the model was placed accurately in the direction of the flow, 
the static pressure on the surface of the body 0:2 in. upstream of the base was 
measured. Experiments showed that conditions at this distance from the base 
were not influenced by the expansion fan at the corner of the model. It was found 
that the upstream influence of the expansion extended for about 0-043 in., i.e. 1:28,, 
where 5, is the boundary layer displacement thickness. Measurements around the 
body were taken by rotating the end piece of the body. They showed that the 
Mach numbers above and below the model were about 0-01 higher than the Mach 
number on either side. 


Static pressure measurements were made to determine the Mach number 
distribution along the centre line of one tunnel side wall and in a vertical plane 
through the axis of the model along a line 0:2 in. upstream of the leading Mach 
line of the expansion fan. The distribution obtained along the tunnel side wall in 
the vicinity of the base of the model is shown in Fig. 4, together with the distribu- 
tion in the empty tunnel. The abscissa gives the distance downstream of the base. 
The fall in Mach number just downstream of the base is caused by a weak shock 
wave which is an effect of the model on the characteristics of the flow. It is visible 
near the base of the model in Figs. 5-8. Because the tunnel and model boundary 
layers are very thick, the expansion fan from the corner of the model affects the 
wall static pressure approximately one inch downstream of the base. The variations 
of the Mach number in the upper and lower parts of the tunnel were found to be 
of the same order as the variations along the tunnel. 
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Fig. 4. 


Effect of the model on the distribution of Mach number along the tunnel side wall. 


From the preliminary experiments it was concluded that the model had only 
a small effect on the uniformity of the flow. The introduction of the model 
increased the mean value of the tunnel Mach number from 1-969 to 1-994. 


The boundary layer just upstream of the base of the model was traversed with 
a small pitot tube. It was found that the boundary layer thickness 5 was 0°18 in. 
and that the boundary layer displacement thickness 5, was 0-036 in. 


The Reynolds number based on displacement thickness was 0:9 x 10*. 


5. Effect of Stings on Base Pressure 


Schlieren photographs were taken of the flow at the base of the model. Typical 
photographs which show most of the features pointed out in Fig. 1 are given in 
Figs. 5-8 (see pp. 227, 228). Because of the highly turbulent nature of the annular 
mixing region and the short duration of the photographic exposure, the “ dead air ” 
region is not visible. The progressive reduction in geometrical scale of the “ dead 
air” and mixing regions with increasing sting diameter causes an upstream move- 
ment of the recompression shock wave and a decrease in spread of the recompression 
fan at the base of the shock wave. In Fig. 9 the base pressure measurements are 
shown in coefficient form against d/D. Since the base pressure for d/D=0 is the 
mean of a very large-number of measurements, the curve has been drawn through 
this point. In the majority of cases the base pressure was measured on both sides 
of each sting, both results being shown in Fig. 9. 


The experiments were later repeated with the boundary layer thickened 
artificially. Air from the atmosphere was bled into the model boundary layer near 
the nozzle throat. This increased the boundary layer displacement thickness just 
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Effect of the ratio of sting to base diameter on base pressure coefficient. 


upstream of the base to 0-045 in. and the Reynolds number based on the displace- 
ment thickness to 1:1 x10‘, without appreciably altering the velocity profile. 
There was no significant difference between the base pressures measured in the two 
sets of experiments, but the experimental scatter was larger in the second set. This 
was possibly because of the moist air in the boundary layer. 


6. Discussions of Results 


The curve in Fig. 9 shows the form to be expected from the discussion in the 
introduction. It shows a gradual decrease in pressure coefficient until d/D equals 
0:8. The curve then begins to turn, reaching a minimum when d/D equals 0°85. 
The curve then rises rapidly to zero as d/D increases to unity. 


The first part of the curve corresponds to the fall in base pressure towards the 
two-dimensional value. This part of the curve is extrapolated in Fig. 9 to the two- 
dimensional value’ which would be reached if (D—d)/D could tend to zero with- 
out becoming small compared with 6/(D—d). When d/D equals 0°85, the curve 
reaches a minimum value of C, at —0:186. This is approximately 0-8 of the two- 
dimensional value under similar conditions. The curve then rises rapidly to zero 
as d/D approaches unity and 6/(D—d) becomes large compared with (D—d)/D. 
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Fig. 5. 
Flow at the base of the model, d/D=0. 


Fig. 6. 
Flow at the base of the model, d/ D = 0-039. 
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Fig. 7. 
Flow at the base of the model, d/ D=0°325. - 
Fig. 8. } 
Flow at the base of the model, d/ D -=0-875. 
228 The Aeronautical Quarterly 


* 
“ 
we) 
‘ é 
it 
iy 


EFFECT OF STING SUPPORTS 


Ratios of sting to base diameter as high as 0-4 are sometimes used in model 
experiments. Fig. 9 shows that, at the Reynolds number and Mach number of the 
present experiments, this introduces an error of approximately ten per cent. in 
the base pressure coefficient of an axially symmetric body. Even with the more 
modest value of 0-2 for the ratio of sting to base diameter, the error in the base 
pressure coefficient may still be as large as three per cent. 
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Photoelastic Investigation in Connection with 
the Fatigue Strength of Bolted Joints 


H. T. JESSOP, C. SNELL and G. S. HOLISTER 


(University College, London) 


SumMmMaryY: The elastic stress distribution around a circular hole in a flat bar 
under simple tension, when the hole is filled by a push-fit pin, was investigated 
te photoelastically. Over a range of values of the ratio of hole diameter to 
" width of bar, and for pins of differing Young’s moduli, the effect of the pin ' 

was sensibly the same, namely, to reduce the maximum tension on the hole 
bi boundary by about 15 per cent. as compared with that in the unfilled hole. 


The maximum shear stress on the boundary was unaltered. 


1. Introduction 


The investigation described here forms part of a research programme initiated 
by the Royal Aeronautical Society and outlined in Ref. 1. In this research it is 
proposed to compare the results of fatigue tests on light alloy specimens with the 
elastic stress distribution as determined by measurements on photoelastic models 
of similar form. The first stage of the programme dealt with the effect of hole 
size upon the stress distribution in a rectangular tension bar with a central circular 
hole. The results of the photoelastic investigation of this were published in Ref. 2. j 
The second stage deals with a similar tension bar with circular hole, the hole in } 
this case being filled by a pin with a push fit. The photoelastic results for this stage 
are given here. 


NOTATION 
E  Young’s modulus 
P.q principal stresses 
t, mean tensile stress in undrilled bar 
t, mean tensile stress across section through hole centre 
tm maximum tensile stress on hole boundary 
The stress concentration factor (S.C.F.) is the ratio t,,/t,’, 


maximum tensile stress 
mean tensile stress across section through hole centre ° 
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2. Method 


The photoelastic materials used were Araldite B for measurement of the stress 
differences and Perspex for observation of the directions of principal stresses. The 
test bars were 1 in. wide and 4 in. thick, and the four sizes of hole investigated 
were 4 in., 3 in., } in. and } in. diameter. 


The three materials used for the pins were: — 


E E (pin)/E (bar) 
Araldite 4:5 x 10° Ib./in 1 
Filled Bakelite 9-7 x 10° 2:15 
Brass 150 x 10° 30 


A large number of different materials were tested in an attempt to obtain a ratio of 
3: 1, but the filled Bakelite gave the nearest value to this that could be found. 


The obtaining of correct fit of the pins presented great difficulty. Turning the 
pins to measured diameters proved inadequate, and it was found necessary to 
check the fit in each case by “feel.” The criterion eventually adopted was that the 
pin should 


(i) push and turn easily in the hole; 
(ii) not fall out when lightly tapped; 
(iii) cause no optical stress effect in the bar when inserted. 


It was appreciated that with pins varying in diameter from } in. to 4 in. this would 
not necessarily indicate “equivalent” degrees of fit, but it was considered the 
only condition which could be attained and repeated with a reasonable degree of 
accuracy. For the } in. hole pin an increase in diameter of the order of 0-00015 in. 
produced a “tight” fit with appreciable stress round the hole. 


3. General Results 


3.1. EFFECTS OF VARYING APPLIED LOAD 


In the original scheme, it was proposed to take observations for two loads which 
produced “lift” of the hole from the pin and two which did not produce lift, for 
each of the three plugs. 


In this push-fit case it was found 


(i) that lift appeared to occur at very small values of load, and 


(ii) that within the limits of error of the observations the elastic stress 
distribution was not altered by alteration of load. 


3.2. EFFECTS OF VARYING THE E OF THE PIN 


It was also found that the stresses were the same with the three different pins 
for all loads up to the safe maximum for the model. 
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Fig. 1. 
Photograph of “ lift” in 4 in. diameter hole filled by a push-fit Bakelite pin. 
Load on bar 110 Ib. 


It would appear that the strains produced by the small transverse pressures on 
the pin were not large enough to allow differentiation between the three cases. 


The results which follow, therefore, apply to all four hole sizes and to the 
three different pins. 


3.3. OBSERVATIONS OF “LIFT” 


The extent of lift occurring was estimated by direct observation by eye, and 
by observation of an enlarged projected image (magnification 7 : 1). 


Figure 1 shows a photogfaph of the lift in the case of a 4 in. diameter hole 
under a load of 110 lb. Fig. 2 shows the angle of lift observed under different loads. 


In view of the scatter in these values and of the difficulty of detecting the extent 
of the lift, especially under small loads, the variation with load is not considered to 
have any great significance. 


It will be noted also that the value of lift obtained for any given load in the 
case of the brass pin lies between the corresponding values of lift for the Araldite 
and Bakelite pins. The results shown are for the 4 in. hole in which the effect was 
most easily observed. 
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Fig. 2. 


Variation in observed angle of “lift” with load. 


It is thus concluded that the observed difference between the three cases may 
well be attributed to small differences in fit between the respective pins and to the 
nature of their finished surfaces. 


4. Detailed Results 


4.1. THE STRESSES ON THE TRANSVERSE SECTION THROUGH THE CENTRE OF 
THE HOLE 


Figures 3 to 6 show the values of the stress differences and of the separate 
stresses for the four hole sizes, in each case under a load which produces two 
fringes relative retardation in the undrilled bar. Values of the stress differences 
for the unfilled holes have been included for comparison. 


It should be noted that in obtaining the graphs of the q stresses it was 
impracticable to carry out the integration over the range aa’ (Fig. 3). Integrations 
were carried out from A to a, and from A’ to a. The resulting curves were then 
smoothly connected so as to satisfy the condition that q had a maximum value at 
the cupic point C whose position was known. The value of the q stress at the 
edge of the hole is therefore not very well determined, and may be subject to an 
error of the order of 0-25 wavelength. 
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4.2. THE STRESSES ON THE VERTICAL SECTION THROUGH THE CENTRE OF 
THE HOLE 


Figures 7 to 10 show the corresponding stresses at points on the vertical axis 
through the hole, in the filled hole case. 


The corresponding values for the unfilled hole do not differ materially from 
these. 


4.3. STRESS CONCENTRATION FACTORS 


Figure 11 shows the stress concentration factors with respect to hole diameter 
for the horizontal section. 


The individual readings are marked by a vertical line, the length of this line 
being an estimate of the probable error incurred in determining the individual 
values. It does not, however, take into account the errors due to differences of 
fit between the respective pins. 


5. General Conclusions 


The insertion of a push-fit pin into the hole under investigation for the case 
where the bar alone is subject to tension would appear to have the following 
results : — 


(i) The peak shear stresses on the hole boundary showed no significant 
change from those in the unfilled hole tests. 


(ii) The peak tensile stresses on the hole boundary were in all cases reduced 
by the presence of the push-fit pin. 


(iii) Any effect due to the elastic modulus of the pin is negligible in com- 
parison with the stress differences produced by very small variations in 
the fit of the pin. 


(iv) The angle of lift does not appear to vary to a very great extent with 
variations of load and the stress distributions are not sensibly affected by 
these variations of load. 


(It should be noted that the stresses in the uniform section of one of the Araldite 
bars were of the order of 960 Ib./in.?, and this would produce a strain of 2 x 10-* 
in the uniform portion. 


The stress to produce equivalent strain in a light alloy bar would be of the 
order of 20,000 Ib./in.2.. Thus the strains in the photoelastic models covered a 
wider range than would normally be encountered in the prototype. The effects due 
to the push-fit pins in the prototype would not therefore be expected to be greater 
than those in the photoelastic model.) 
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Ye 4 
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Fig. 11. 


Stress concentration factors for filled and unfilled holes. S.C.F.=1,,/1,’ where t,,=maximum 
tensile stress on hole boundary, t,’=mean tensile stress across section through hole centre. 
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REVIEW 


Review 


“Linearized Theory of Steady High-Speed Flow,” G. N. Ward. Cambridge University 
Press, 1955. 243 pp. 30s. 


Extensive use has been made of the linearised equations of motion in the practical 
development of high-speed (and especially supersonic) aerodynamics. In many cases it 
is only by using these linearised equations that any solutions can be obtained, and the 
limitations introduced by the approximations made in the theory must be accepted. 


Professor Ward’s book will be welcomed by mathematicians and theoretical aero- 
dynamicists working in this field. The book is not intended to be a comprehensive 
treatise, covering all aspects of the subject, but it does give a logical development of 
the basic theory, with a careful discussion of the assumptions and approximations, 
and then goes on to consider some methods of solving particular problems. 


The first part of the book deals with the general development of the theory, 
considering the mathematical implications in some detail, and including an account of 
flow-reversal theorems. The second part starts with a short discussion of subsonic flow, 
and then goes on to consider supersonic flow past nearly plane wings and the theory 
of conical fields in supersonic flow. The last chapter of this part describe§ the use of 
operational methods in problems of supersonic flow. Finally, in part three, there is 
an account of slender-body theory. 


The author has made extensive use of vectorial notation in developing the theory, 
and he draws attention to the advantages of this in revealing the physical meaning of 
the equations. Since the book is essentially one for mathematicians, and not for 
engineers of only modest mathematical ability, the use of vectorial notation should not 
cause any difficulty among serious users of the book. If any further warning is needed 
that the book is for mathematicians only, it may be sufficient to remark that it is a 
slightly revised and shortened version of an essay submitted in competition for the 
Adams Prize for 1949-50 in the University of Cambridge. The author remarks in the 
preface that this explains the “ strong personal bias ” in the selection of material. Since 
Professor Ward has himself made important and extensive contributions to linearised 
theory, this “ personal bias” will be welcomed by his readers. 

W. A. MAIR. 


Addendum — May 1955 


After the printing of the paper by E. W. Graham, B. J. Beane and R. M. Licher, 
“The Drag of Non-Planar Thickness Distributions in Supersonic Flow” in the May 
1955 issue of The Aeronautical Quarterly (Vol. VI, p. 99) the authors found that the 
idea of the central body plus cylindrical shell with zero drag, attributed in the paper to 
Ferri (Ref. 5), was in fact proposed earlier by Ferrari. The reference to this earlier 
work is:— 

FERRARI, C. Campi di Corrente Ipersonora Attorno a Solidi di Rivoluzione, 

L’Aerotecnica, Vol. XVII, No. 6, pp. 507-518, 1937. English translation: Super- 

sonic Flow Fields about Bodies of Revolution. Brown University, Graduate 

Division of Applied Mathematics, Translation No. A9-T-29, prepared for Analysis 

Division, Intelligence Department, Air Materiel Command, Wright Field, Dayton, 

Ohio, 1948. 
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Note on a Null Method for the Direct 
Measurement of Pressure Coefficients 


L. H. TANNER, B.A. 


(National Physical Laboratory) 


SuMMaRY: Values of the pressure coefficient C, between + 1 and - 1 can be 
measured directly in incompressible flow by comparison with the known 
pressures on a circular cylinder placed in the tunnel working section. The 
development of the method and the apparatus required are described. 
The method gives more consistent and accurate readings than can easily be 
obtained with normal manometer measurements, particularly when the 
pressure differences are small and the tunnel speed not very steady. Since 
the coefficients depend on Reynolds number only, less precise measurement 
and control of tunnel speed are required, unless the coefficients to be 
measured are subject to critical scale effects. Possible ways of extending 
the range of coefficients which can be measured are suggested. 


1. Introduction 


In almost all aerodynamic experiments which involve pressure measurements, 
a pressure ratio, or a ratio of pressure differences, is required, rather than the actual 
value of the pressure. In high-speed work the quantity required is sometimes the 
ratio of the pressure on the model to the stagnation pressure, but in low-speed work 
pressure coefficients—ratios of two pressure differences—are almost always used. 


In view of this it seems surprising that methods for measuring these ratios 
directly are so rarely used. The only instruments for this purpose in general use 
seem to be the various types of Machmeter; these are normally used only for setting 
or measuring the free stream Mach numbers of wind tunnels or aeroplanes, not for 
routine pressure measurements. The reason why such pressure ratio instruments are 
so little used clearly lies in the simplicity and convenience of liquid manometers, 
particularly of multiple manometers for rapid reading of pressures at many different 
points. Such manometers have, however, some disadvantages which are worth 
remembering. With an ordinary multi-tube manometer, each result involves three 
meniscus readings, each of which contains an error. This.error arises from the 
reading error and from the effects of surface tension, clamping the tubes, before 
reading, and errors in the tube position in the case of an inclined manometer. The 
accuracy depends largely on the magnitude of the pressure differences to be 
measured and becomes poor if these are small. 
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is only by using these linearised equations that any solutions can be obtained, and the 
limitations introduced by the approximations made in the theory must be accepted. 


Professor Ward’s book will be welcomed by mathematicians and theoretical aero- 
dynamicists working in this field. The book is not intended to be a comprehensive 
treatise, covering all aspects of the subject, but it does give a logical development of 
the basic theory, with a careful discussion of the assumptions and approximations, 
and then goes on to consider some methods of solving particular problems. 


The first part of the book deals with the general development of the theory, 
considering the mathematical implications in some detail, and including an account of 
flow-reversal theorems. The second part starts with a short discussion of subsonic flow, 
and then goes on to consider supersonic flow past nearly plane wings and the theory 
of conical fields in supersonic flow. The last chapter of this part describés the use of 
operational methods in problems of supersonic flow. Finally, in part three, there is 
an account of slender-body theory. 


The author has made extensive use of vectorial notation in developing the theory, 
and he draws attention to the advantages of this in revealing the physical meaning of 
the equations. Since the book is essentially one for mathematicians, and not for 
engineers of only modest mathematical ability, the use of vectorial notation should not 
cause any difficulty among serious users of the book. If any further warning is needed 
that the book is for mathematicians only, it may be sufficient to remark that it is a 
slightly revised and shortened version of an essay submitted in competition for the 
Adams Prize for 1949-50 in the University of Cambridge. The author remarks in the 
preface that this explains the “ strong personal bias” in the selection of material. Since 
Professor Ward has himself made important and extensive contributions to linearised 
theory, this “ personal bias” will be welcomed by his readers. 
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